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ABSTRACT 


A method for passive ranging based on the principle of triangulation is considered 
In the basic triangulation scheme that is a single baseline model the precision in the bearing 
readings can be related to the precision in the range estimation For some target orienta¬ 
tions the precision in the triangulated target range is completely lost This phenomenon is 
known as "geometric dilution." A proposed orthogonal dual baseline scheme eliminates 
the "geometric dilution" effect. The performance of each of the two orthogonal baselines 
depends on target orientation. For specific target orientations the triangulation range mea¬ 
surements for the two baselines are equivalent. The dual baseline scheme would require 
"smart electronics" which would switch between baselines at crossover points in the range 
estimation precision It is shown that the crossover points depend primarily on the ratio of 
the two baselines. A general expression for the maximum triangulation range consistent 
with limitations in minimum tolerance precision in range estimation is derived. The depen¬ 
dency between maximum range and target orientation are presented in polar form. Limita¬ 
tions in the dual baseline model due to physical limitation created by the optical horizon 
are also considered. 
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I. INTRODUCTION 


A. REVIEW OF PRIOR WORK 

Recent technological developments can turn an active radar system into a target for 
anti-radiation missiles (ARM). In essence, a system designed to provide security against 
missile attacks can become a security hazard. This thesis will present a methodology of 
passive ranging based on triangulation and the use of IR sensors which could be imple¬ 
mented on various kinds of military platforms. IR sensors provide for passive surveillance 
and target identification and tracking, but targeting and fire control require a minimum 
range to target 

Passive ranging is the Holy Grail of the defense community. The passive interroga¬ 
tion of a scene greatly restricts the domain of applicable countermeasures. Therefore, with 
passive ranging the security of the intelligence operation is less likely to be compromised. 

Various passive ranging schemes based on the radiating characteristics of a target 
have been proposed These can be classified under a variety of descriptors such as time 
difference of arrival (TDOA), frequency difference of arrival (FDOA), and angle differ¬ 
ence location (ADL) [Ref 1], [Ref 2], The first two schemes, TDOA and FDOA require 
that RF modulation be present on the target's electromagnetic spectrum [Ref 3] Since 
the thrust of this thesis is to focus on a potential ADL technique triangulation, which re¬ 
quires only a target's unmodulated IR signature, no further discussion on TDOA or FDOA 





methods will be included. Several methods for using atmospheric attenuation to predict 
target range passively are covered in the next paragraph 

In the early 1960's several patents were approved which related IR signal's atmos¬ 
pheric attenuation to range. The first of these was based on the atmospheric absorption of 
radiation from targets which have strong water vapor or carbon dioxide bands [Ref 4] 
The scheme includes signal chopping, optical beam splitting, and controlled absorption of 
the CO, band in a reference path A nulling circuit, which uses a partially transmitting 
wedge driven by a servomotor, balances the two signals. The error signal is then directly 
related to the amount of energy in the CO^ band of the original IR signal received This er¬ 
ror signal can be related to the length of the attenuation path, i .e., range. A second patent 
[Ref 5] is again based on the principle of preferential atmospheric absorption of IR radia¬ 
tion for some bands. However, in this case the authors proposes that if a target’s total and 
filtered radiation intensities are known, and if the effective radiation coefficients for total 
and filtered bands could be ascertained, then ratios of measured intensities can be related 
through analysis to the unknown range The inventors of this patent also use a chopper 
and comparison circuitry. Recent NFS theses have examined the passive ranging problem 
using similar concepts [Ref 6], [Ref 7], [Ref 8]. 

The principle of the parallax effect, as applied in the optical range finder [Ref 9], has 
been modified to process IR signals. In this scheme, introduced as a patent [Ref 10], the 
measurement of a parallax angle between two imaging sensors and a known baseline will 
permit calculation of the range to a target. Another 60's patent [Ref 11] suggests a model 
for calculation closure time, the approximate time it will take one body to collide with 
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another. This parameter depends on range. The inventors show that the closure time can 
be calculated from the time rate of change of the target's passively detected radiant energN' 
They provide a circuit for this application. 

With the development of IRST passive sensors [Ref 12], renewed interest in .ADL 
methods has occurred. In particular these devices are available on a rotating (30 rpm) plat¬ 
form and can predict bearing azimuth data on low flying targets A platform for an IRST 
device is shown in Figure 1 The bearing resolution of such devices is classified; however 
an interesting mathematical model for predicting resolution for targets in the presence of 
noise has been proposed [Ref 13]. 



Figure 1: ERST sensor platform (After Ref 12) 








Most of the ADL methods are based on using two sensors, i e , a single baseline ap¬ 
proach (SBL) [Ref 14], [Ref 15] In the late 70's a multiple sensor scheme based on a te- 
chique known as circulation [Ref 16], which employs an interactive least squares 
algorithm, has been investigated Although the method eliminates the problem of "geomet¬ 
ric dilution" the scheme can exhibit convergence problems. More recently a preliminary 
study on an analytic dual ba.seline method has been introduced [Ref 17] This description 
is based on an approximate model for estimating the range sensitivity Further develop¬ 
ments on this dual baseline model will be presented in this thesis 

B. THESIS OVERVIEW 

Chapter I covers the literature review and the thesis overview. In Chapter II the sin¬ 
gle baseline analysis is developed, and the lack of precision in a triangulated range is re¬ 
lated to the lack < f precision in the bearing measurements. It is shown that the precision of 
the triangulation method for targets alligned with the single baseline is completely lost 
This problem is known as "geometric dilution " Chapter III explores the applicability of a 
dual baseline scheme in order to circumvent this problem. In Chapter IV the dual baseline 
model is applied to characterize a maximum range for triangulation consistent with a mini- 
am precision in the measurement. The development in this chapter leads to a physical 
polar representation for in terms of target orientation Two paradigms, quadrature 
addition and the more conservative sum of absolutes, are considered In Chapter V the 
limitations in range estimation due to earth's curvature are evaluated and included with the 
dual baseline model Some interesting qualifications for predicting the characteristics of 
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the polar plots studied in Chapter IV are provided in Chapter VI Chapter VII provides 
the conclusions 


Appendix A is the list of symbols provided for convenience in alphabetical order In 
order not to complicate the main text some of the more mathematical derivations have 
been relegated to the Appendices B, C, and D In Appendix B details on approximations 
made on the single baseline model are presented. In Appendix C, an exact mathematical 
expression or the single baseline triangulated range in terms of IRST sensors bearing is 
derived. In Appendix D two paradigms for error estimation are analyzed. It is shown that 
the sum of absolutes estimation is more conservative than estimation via quadrature addi¬ 
tion Appendix E is a collection of the programs used in generating the plots seen in the 
thesis 
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11. BACKGROUND ANALYSIS 

A. THE TRIANGULATION CONCEPT 

The eye-brain system undoubtedly works in a fairly sophisticated manner in order to 
extract information from parallax. The geometric considerations of the corresponding pas¬ 
sive sensor problem are fairly simple. In principle, the range prediction can be made with 
100% accuracy; in practice the limitations in the accuracy of the measurements in orienta¬ 
tion bearings introduce an error in the range prediction. The sensitivity for the range pre¬ 
diction to variation in the bearing accuracy is dependent on target orientation, the distance 
between the points of observation or baseline, and the range to the target Currently IRST 
bearings do not include elevation data. For targets which appear coming over the horizon 
this should not represent a significant deficiency in the model. Therefore, the model to be 
presented predicts horizontal range, not the slant range. 

B. SINGLE BASELINE MODEL 

A single baseline model, involving two sensors, located at the ends of a baseline, can 
be analyzed for conditions sufficient foi triangulation. Starting from Figure 2 the points at 
which bearing measurements are made, A and C, are located at a distance R, and IL re¬ 
spectively from the target at T. Although no measurements are made at point B it serves 
as a convenient symmetrical reference between observation points A and C. The range R 
is defined as the distance from point B to the target. The characteristics of the single 
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baseline model in terms of symmetric variables R and 0 greatly facilitates the dual baseline 
analysis presented later. The distance between the points A and C is the baseline (BL). 
The initial objective of this analysis will be to relate the range R to the measurables in the 
problem, i.e., baseline (BL) and orientation (0, and 02). The next step in the procedure will 
be to perform a sensitivity analysis with respect to variation in the orientation variables, 0, 
and 02 - 

The significant assumptions of the model are; 

(1) There is no error in the baseline measurement; 

(2) The tolerances in the bearing measurements are known, and will be denoted 
A0, and A02; 

(3) The discrepancy between horizontal path range and slant-path range is not 
significant. 

1. A single baseline analysis 
From the Figure 2, the apex angles satisfy the following relation 

03 = 02 - 0, = 0, + 0, (1) 

Another relation between the apex angles and 0,, ©2 is obtained from the Law of Sines 

sin 04 = ^ sin 02 (2a) 

sin 05 = ^sin0| (2b) 
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Through expansion of (1) it follows that 


sin (62 -6i) = sin64Cos65 + cos64sin0.<; 


(3a) 


and 


sin ( 62 - 61 ) = ^ 


sin02jl -^^sinOi j +sinei - ^^sinS- 


(3b) 


Solving (3) for the range R, it follows that 


sin 02 ^ 

fl-^^sin 0 ij +sin 0 )^ 

(l-(isin 02 ]‘ 

sin (02 - 0 )) 


(4) 


An exact expression for R/BL(0,, 0^) can be obtained using symbolic mathematical proc¬ 
essing. See Appendix C. The immediate objective is to obtain an expression for the range 
sensitivity in terms of the bearing angles 0 , and 0 ^, which will be in the form 


^ -/i(0|,02)A01 +/2(01,02)A02 


(5) 


In order to not distract the reader, the mathematical details in obtaining the range sensitiv¬ 
ity appear in Appendix B, and the final expression for it is 
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where it is understood that AO,, and A6, are assumed to be zero mean random variables 
A statistical model for eliminating "systematic error" has been proposed using a scheme 
based on a method known as circulation [Ref 16], Furthermore, the ratio "a" of rms val¬ 
ues for A0 t and A0, can be used and the final expression for the normalized variation in R 
can be written as follows 


I-(|j^sin02 j (cos0|+«in9icol(62“0i)>+5m8j^l-^|^$in0i j cot(02-0i>— 
Sin 02 j +sm 0 i il-[|j^sin 0 j) 


nneicme. 


l-[|ysin 0 ij {cos02-sin03cot(0j-0i))-«ui0i /l-[§sin02j col(02-0|)+-^^^- 

_ 

nn02 l-[||^iin 0 ij ■MU10IJl^^Jin02j 


• X2 • O 


where %] and X 2 are independent random variables having zero mean and unity standard 
deviation, and 



(8) 
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It is noteworthy that (7) depends on R, 6,, and 6, instead of the symmetrically defined 
variables R and 6. Additional relations are needed at this point in order to express (7) and 
(8) in terms of the symmetrically defined variables R and 6. From the Law of Cosines, 
applied to Figure 2, it follows that 


"i )' = (i) +(0 




From the Law of Sines, applied also to Figure 2, it follows that 


■ n • n(«) 

Sin 02 = sinOjT^ 


(10a) 


sin6i =sin6“rT 


(10b) 


and, again from the Law of Sines applied to the same figure 


sin(02 -0i) =sin(04 +05) = |^sin0i = |^sin02 

K\ «2 
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Finally, along with the trigonometric rule 


COtX = ' 


substitution of (8), (9), and (10) into (7), which is the normalized sensitivity of the range, 
limits the dependence of variables to R and 0 only. Although this expression is not shown 
explicitly in this form, it is readily evaluated in its implicit form via the substitution 
procedure. 

2. Approximate single baseline analysis (R»BL) 

To appreciate the factors that primarily limit the accuracy in the predicted 
value for R it can be instructive to consider a simpler formulation of the problem. Consis¬ 
tent with the conditions of assumption 3 the ratio of the range to the baseline can be con¬ 
sidered much bigger than one (R/BL»1). After application of this assumption to (3) the 
second order terms in BL/2R can be neglected, leading from (4) to 


p __ sm j ^ 

^- 2 ]^ sm{t»2-0,) ) ~ 2 


and the normalized variation in range would be 


Sin 02-^0101 


■cot(e,-e,))(A9,) + (:;K5^ -cot(e:-e,)j(Ae2) 




Equation (6) representing the normalized variation in R can then be expressed 


Ag 

K AH 


sinH:-KinB 


+ cot (02 - e,) ] • Xi + - cot (6: - e,)) • X: • 


(15) 


where again x, and X: are independent random variables having zero mean and unity stan¬ 
dard deviation, and "a" is defined by (8). Equation (15) can re-expressed in terms of R 
and 0 after taking into account (9) through (12). Again the details appear in Appendix B 
3. Approximate single baseline analysis (R»BL, R=R,= R^) 

From the Law of Sines, Equation (8), it follows that 

BL 

^ 1 ~ ,/> , sin 01 /1 z \ 

sin(e2-«i) ( 16 a) 


and 


^2=-ifti;ysin02 (16b) 

Working with (12), and (13) the reasonable assumption R=R, =R 2 can be made and (14) is 
easily converted into the approximate form 

^ = ^^cot0,+cot(02-e,)j(Ae,)+(^|cote2-cot(e2-0,)j(Ae:) (i7) 
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after noting that R = R\ =R 2 implies 


sui6| = sin6: HsinO 


( 18 ) 


Furthermore from (16) 


sin (02 -0i) = ^sin6 


(19) 


and substitution of (18), and (19) into (17) leads to 


R 



l-(^sin8) 

f 


1 , 1 

(Ae,)+ 

1 ^l-sin^0 ^ 

2 sme 

fsme 

2 sm0 


V 

J 



/ 


(Ae,) 


( 20 ) 


Using the ratio "a" of A02 and A0, applied to (19), it follows that 



f 

(77 7= > 



{ - 

1 Jl-**"'* tl 


•Xi + 


R A9 J 

2 sin0 

&su,9 

2 sin9 



J 


k / 


( 21 ) 


This expression represents the normalized sensitivity in terms of the correct symmetrically 
defined variables, R and 0. 

4. Comparison of single baseline models 

In order to continue the triangulation analysis a comparison between the three 
single baseline models. Equations (7), (15), and (21), has to be made Due to very com¬ 
plicated expressions that stand for the three models a reasonable comparison of them can 
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complicated expressions that stand for the three models a reasonable comparison of them 
can be made through graphical representation. In order to compare the three approximate 
forms for the range sensitivity, the independence of the random variables Xi and X- should 
be considered. Since all three approximate expressions for the range sensitivity are of the 
form 


TX? = Uerml) • x, + (term!) ■ X: 


( 22 ) 


the natural method of comparison is to take 

I I = 7(/enwI)^T7/crw2^ (23) 

or equivalently the coefficients of the sensitivity add in quadrature. This, as suggested in 
Appendix D, may not be the most conservative estimate for the sensitivity However for 
purposes of evaluating the accuracy of the approximate model (21) relative to (7) and (15) 
it is adequate. The three approximation levels of the range sensitivity expressed in Equa¬ 
tions (7), (15), and (21) are compared in the Figures 3, 4, and 5 for three values of the 
range to baseline ratio, and for ratio of bearing errors 0, 1 and 2. 

5. Observations 

Scrutinization of Figures 3, 4, and 5 indicates only small discrepancies between 
the approximate solutions. It is concluded that the approximations of Equation (21), 
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characteristic of all the single baseline figures are worth noting. First the range sensitivity 
is highest for angles close to 0" and 180“. This is the well known "geometric dilution" [Ref 
14] that is characteristic of all the single baseline methods. In terms of measurements a 
high normalized sensitivity would result in a high level of inaccuracy. This provides the 
motivation for exploring the dual baseline method which is the focus of the later chapters 
Because the range sensitivity normalization is done with the respect to AS, the effect of in¬ 
creasing "a" is to amplify the term proportional to AG,. Therefore, as seen from the fig¬ 
ures, there is an overall increase in the normalized range sensitivity. Lastly, the farther the 
triangulated target is from the baseline the more inaccurate the measurement. This is also 
confirmed by the increase in the normalized range sensitivity of Figures 5a, 5b, and 5c 
over the respective counterparts of Figures 3, and 4. 
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T 



Figure 2: Geometry for the single baseline model 


16 





Idella R / (R * delta lhela)( 



Theta in degrees 


Figure 3a; Single baseline model sensitivity curves for R/BL=10 & a=0 
(via program in Appendix E. A) 
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comparison of single baseline sensitivity curves 
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o for R»BL approximation case (equation 15) 

+ for R»BL & R-R1-R2 approximation case (equation 21) 
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Figure 3b : Single baseline model sensitivity curves for Rj'BL=10 & a=l 
(via program in Appendix E. A) 
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Figure 3c Single baseline model sensitivity curves for R/BL=10 & a=2 
( via program in Appendix E A) 
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comparison of single baseline model sensitivity curves 
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0 for R»BL approximation case (equation 15) 

+ for R»BL & R-R1 “R2 approximation case (equation 21) 
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Figure 4a . Single baseline model sensitivity curves for R/BL=100 & 0=0 
(via program in Appendix E. A) 
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Figure 4b Single baseline model sensitivity curves for R>'BL= 100 & 1 
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Figure 4c : Single baseline model sensitivity curves for R/BL= 100 & a= 
(via program in Appendix E A) 
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Figure 5a Single baseline model sensitivity curves for R/BL= 1000 & ci=0 
(via program in Appendix E. A) 
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Figure 5b ; Single baseline model sensitivity curves for R/BL=1000 & a -1 
(via program in Appendix E. A) 
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Figure 5c : Single baseline model sensitivity curves for R/BL=1000 & a=2 
(via program in Appendix E. A) 
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m. DUAL BASELINE MODEL 

A. DUAL BASELINE MODEL ANALYSIS 

Figure 6 shows a dual baseline model. The two baselines of the model are taken to 
be orthogonal and they are referred to as the long base line (LBL) and the short baseline 
(SBL) system. The features and the relations of the single baseline model analysis are ap¬ 
plicable here but the superscripts S and L will be used to distinguish the two systems. 
Furthermore, some new relations between the angles are needed, in order to express the 
short baseline (SBL) model in terms of the same variables R and 6, as in the long baseline 
(LBL) model. 

B. GENERAL CONSIDERATIONS 
From the Figure 6 it is clear that 

cos 0^ = cos (6 - 90’) = sin 0 ( 24 ) 


similarly 


sin0^ = -cos0 


Therefore from Equations (9), (24), and (25) 

[m] =[m] +(i)+(2lz:)cos0 


(25) 


(26a) 
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(26b) 




Similarly, from Equations (10), (24), and (25) it follows that 


sin of = sin 0^^^ 


LBL 



sin0| = “COS0 


(sir) 


sinSf =-cos 07 ^ 


and from Equations (11), (24), and (25) 


sin ^02 - 0t 


] 


_ LBL nL _ LBL ■ nL 

-■^sm0, = —sin02 


sin 



_ SBL • ns SBL ■ n 

= -^sin0| =-^sm0 




(26c) 


(26d) 


(27a, b) 


(27c, d) 


(28a) 


(28b) 
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For computation of j 

(S:)~^(z^] (S)^(m] 01 ^ei 02^02 (29a, b, c, d) 


and follow substitution (26a, b) and (27a, b). For computation of 


f ~ t' 


let 



01—>0‘^ 02 —> 02 (30a, b, c, d) 


and follow substitutions (26c, d) and (27c, d). In Appendix B expressions for j 
and j are shown for the various levels of approximation introduced in the previ¬ 

ous single baseline description. As noted in this chapter the approximations R»BL, and 
are quite good, and the dual baseline results (see (21)) for this case are provided 
here for later reference 
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2 sin 6 
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2 sin 9 

(¥•>"») 


•X2«'-(31a) 
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1 /-... \2 ^ 
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fM±Y - 


Ji-(^^cose) 


] ^i-cos^e ^ 

fi-(-^cose] 

[ R A8 J ~ 

2 (-cose) ■ rm,_^„,e)) 
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2 (— COS0) 
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•X^«''■(31b) 
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where X\ ' Xz ^ X:^ 3re independent random variables having zero mean and unity 

standard deviation, and 


= 


a^ = 


(32a, b) 


C. SENSITIVITY CURVES FOR THE DUAL BASELINE SCHEME 

The expressions^—^ j , and in (31) represent the normalized sensitivity 

for the long and the short baseline system respectively, in terms of the symmetrical de¬ 
pendent variables, R and 0. It is obvious from (31) that the sensitivities of both long and 
short baseline systems depend heavily on the terms (R/LBL), {R/SBL\ AO,'-, Ae,^ AB,'", 
and A0/ To illustrate the dependence of the sensitivities to the above mentioned terms 
some examples will be introduced with the use of (31) which, after (23), becomes 
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(33a) 



1 Vi-cos^e y 

2 (-COS0) I SSL 






I » 


(-cos 


e,) 




J__ 

2 (-cos 0) 




LV 


( SSL 
V K 


(-COS0) 


■ a’ 


(33b) 
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where as previously described the terms are added in quadrature However geometrical 
considerations, described in Appendix D, suggest that a more conservative gauge for 
range sensitivity is based on addition of absolutes 

+ \iterm2)\ ( 34 ) 

under this assumption Equation (31) would be evaluated according to 








( 




+ 


] ^l-sin^0 

(Ae j 


2 sin 6 



2 sin 6 



(35a) 


f AR 1 



-t- 



( R A9 j 


2 Mose) ' (^^(^ose)) 


2 (^osB) (^(^ose)) 


(35b) 


It can be seen from Appendix D that for the discussion given, (35) for r/ = </= 1 is nomi¬ 
nally more conservative than (33). Furthermore for c^= 0 , Equation (33) and Equa¬ 

tion (35) are equivalent. 

D. EXAMPLES 

The examples discussed in this chapter will be based on the sensitivity expression us¬ 
ing quadrature addition. Equation (33). The concepts applied are equally applicable to the 
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more conservative equation (35). MATLAB programs for evaluating (33 

Appendix E As a first example take 

) are provided in 


ir = 600 

(36a, b) 

=2 


(36c, d) 

After substitution of the data into (33) and plotting it the results are shown on Figure 7. 

As a second example take 

^=1000 

Wt = 5000 

(37a, b) 

a^ = 1 

a^= 1 

(37c, d) 


The plots following substitution of (36), and (37) into (33) are shown in Figure 8. In 
both Figures 7 and 8 the performance of the short baseline system will exceed that of the 
long baseline for angles close to 0° and 180”. For Figures 7 and 8 those crossing points are 
approximately +17” and +15", respectively. To appreciate the utility of the analysis take a 
simple example (i.e., that one sensor on each baseline has no bearing error,while the bear¬ 
ing errors of the remaining sensors are equal). Then the situation can be described as 
follows: 

A0j=Aef = Ae (38a) 
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ASj = ASj = 0 


(38b) 


thus, the long and the short parts of (33) after setting = 0 become 


^ 1 


M_L _ 

/f 46 J 2 sine 






(39a) 


M_L 'l I i/<-co^ 

/? 46 j 2 cose 


^ Ji-{^cose) 


-cose 


(39b) 


which permit the plotting of the sensitivities versus 6. This was done for the conditions 

^= 1000 and ^ = 5000 (40a, b) 

and the results are represented on Figure 9. Note, as expected, the LBL does better over 
a wider range of angles. Around 6=170” and 0=10” there is a crossover point for which 
the short baseline system would produce more accurate range measurements at angles 
6<10'’ and 0>17O'’. For illustration of the use of the curves assume that the bearing resolu¬ 
tion is A6‘'| =A0^2 =A6 = 0.1 mr and the range to the target is 20 Km. Note according to 
(39) this implies long and short baselines of 20 m and 4 m, respectively For 0=90°, the 
LBL best case, the uncertainty in the long baseline range measurement from (39a) is 

=/?• A6tanv|/^ (41) 
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where AR^ is the uncertainty in the long baseline range, and 


\|/^ =atan 

Then, the range measured by the long baseline will be 

RLas=R±^ieas (43a) 

R^^=20 Km ±2.0 Km (LBL) (43b) 

For 0=0“ (or 180°), the SBL best case, the uncertainty in the short baseline range measure¬ 
ment from equation (39b) is 

^ieas = ^ • Ae • tan (44) 

where AR® is the uncertainty in the short baseline range, and 




Thus, the range measured by the short baseline will be 


(45) 


Rieas =R±ARLos 


(46a) 


R'me« 20 Km ±]0 Km (SBL) 


(46b) 




Although the inaccuracy in the measured value for R may be tolerable in the LBL case, the 
SBL range estimate is probably not acceptable Figure 9 shows what happens to the nor¬ 
malized range sensitivity versus 6 for fixed range. For the dual baseline scheme the meas¬ 
urement accuracy is at a minimum at the crossover points. For the conditions taken this 
occurs at 6=+10“ and 6=-i-170“. 

E. RANGE LIMITATIONS FOR FIXED TARGET ORIENTATION 

Rang"! sensitivity curves plotted versus range for the principal directions 6 = 0“ and 6 
= 90“ could be applied to predict maximum bearing error and/or minimum baselines in or¬ 
der to accurately triangulate targets at desired distance limits. Under the assumptions in¬ 
troduced by (38) the normalized range sensitivities are given in (39), which for the single 
baseline model becomes 



1 ^l-sin^(8) 

2 sin (9) 




1(0)1 


(fsm(e)) 


(47) 


Equation (47) is plotted with R/'BL=\QQ to R/BL=\000, and for angle 0=90", 45",and 22 5" 
on Figure 10. In order to clarify the application of Figure 10 to dual baseline triangulation, 
an example is presented based on the following conditions 

LBL=20.0 m SBL=4.0 m (48a, b) 

0=90“ A0=O. 1 mr (48c, d) 
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and precision in range < 5%, that means 


f < 0.05 (48e) 

From (48d), and (48e) the range sensitivity is evaluated as 

(fi)<500 (49) 

and therefore from Figure 10 the range to baseline ratio is found to be 

K S 550 (50) 

This condition is applicable to both the long and short baseline In both cases the 90" con¬ 
dition (48c) imposes the interpretation that the target appears along a line which is perpen¬ 
dicular to the baseline. 

For the LBL case it follows from Equation (48a) that 

^ 550x20w = 1 IAtw (51a) 

and similarly for the SBL case 

Rsbl X 4.0 m = 2.2 km (51b) 
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This demonstrates that the target must be significantly closer to be triangulated with 5% 
accuracy if the oncoming direction is perpendicular to the short baseline Nonetheless, as 
previously mentioned, the closest unresolved range will occur for the previously noted 
crossover angle. To demonstrate the relevance of this issue. Figure 11 shows the plot of 
(47) for the angles 0=90”, 0=8O”,and 0=10°. The angle 10° was chosen because it repre¬ 
sents the approximate crossover point for the conditions (38). After referring to Figure 11 
it is seen that the performance for the 0=10° case is about six times worse than the 0=90" 
case. Using the LBL baseline (51a) the range detected by the long baseline would be 
approximately 


/?<1.8A7w ((^ 10“) (52) 

which is, as expected, worse than either of the baseline performances (51). It should be 
noted that for a crossover point the long base line and short baseline performances are the 
same This can be checked by using the 0=80” curve with the SBL data. 

F. COMMENTS ON THE EXAMPLES AND ANALYSIS 

The primary objective of the preceding examples and discussion has been to demon¬ 
strate a feasible modeling scheme for dual baseline triangulation The results of this chap¬ 
ter demonstrate several features of the dual baseline scheme. The problem of "geometrical 
dilution" discussed in Chapter I can be avoided by using two orthogonal baselines In gen¬ 
eral the larger baseline will have a higher measurement accuracy. Because of typical ship 
structure it is assumed that a second baseline would not be conveniently as long The 
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scheme would require some "smart electronics" which would switch between baselines at 
a performance crossover point. As seen from the discussion in section D the maximum tri¬ 
angulation range, for a fixed range sensitivity , can be predicted for particular target ori¬ 
entations, i.e., fixed 0. In the chapter to follow this concept is extended by solving for 
versus 0 under conditions of constant range sensitivity. 
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Figure 8: Dual baseline second example(R/LBL= 1000, R/SBL=5000, a^=l, a"=l) 
(via program in .Appendix E.B) 
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Figure 9: Dual baseline third example(R/LBL=1000. R/SBL=5000, -£^=0) 

(via program in Appendix E. B) 
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Figure 10; Range sensitivity for fixed target onentation (0= 90°, 45°, and 22.5°) 
(via program in Appendix E C) 

42 




delta R/{R‘ di 


range estimation for fixed target orientation 


3000 


2000 



500 600 700 800 900 1000 

range/base-line 


Figure 11: Range sensitivity for fixed target orientation (6= 90“. 80°, and 10“) 
(via program in Appendix E. C) 








IV. RANGE LIMITATIONS FOR FIXED SENSITIVITY 

A. CHAPTER OVERVIEW 

In this chapter a more graphical representation for versus orientation angle 6 
will be developed. For target ranges less than the triangulation estimation will, by as¬ 
sumption, be within a tolerable precision limit dictated by a maximum sensitivity. Two 
paradigms for estimation are considered. The first is based on a quadrature addition of 
terms. In this case an analytic solution can be obtained for the case a=\, i.e., equal preci¬ 
sion bearing measurements. In the second case the addition occurs as a sum of absolutes. 
Here after making a reasonable approximation R/BL »1, a solution for all "a" can be 
obtained 

B. ANALYSIS BASED ON THE SUM OF ABSOLUTES 

Define the normalized range sensitivity for long and short baseline respectively as 





^_L 
R AB 




(53a, b) 


after noting for both LBL and SBL that 


2 iin(8) ^sui(9) 


(53c) 
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and after substituting, (35) simplifies to 




(54a) 




(54b) 


For typical cases R »BL (as previously assumed in Chapter II ) the first terms, propor¬ 
tional to ( 1 - 0 ^^), are very small compared with the second ones, proportional to (l+o^'*), 
so it is possible to neglect them without significant sacrifice in the accuracy. Then (54) be¬ 


comes 


(^$ 111 ( 64 )) 


(1 +a^) = c^ 


(55a) 


(^»ir(9^)) 


(l-l-o^) = c^ 


(55b) 


Solving (55a) for it follows that 


/ \ 2 

[Srj -sin (6)-^ 


LBL ~ t 


(55c, d) 
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Equation (55d) can be used to create a polar plot of the range measured by the long base¬ 
line normalized to the LBL Solving (55b) for SBL, and after (25), it follows that 


f<SBl 

SBL 


cos‘(6 ) 


(i+a-'r 


(55e) 


Bsbl 

LBL 


COS*(0) 





(55f) 


Equation (55f) gives the range measured by the short baseline normalized with respect to 
the LBL for convenient comparison with (54). 

C. ANALYSIS BASED ON QUADRATURE ADDITION OF TERMS 
After the reasonable assumption that a = 1, (33) becomes 



j/' 
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fAff 1 '1^ _ 

1 ^ 




“ sui(9^) 
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2 sin(B^) 
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(56a) 


^_L 

R AB 


J 2 5in(e'') 
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2 sin (9*) 
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(56b) 


After defining 


1 

2 sin(B) 


COt(9) =:^ 


(57) 
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and 
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After referring to Figure 5 and solving (60a) for it is found that 


f<LBL 

LBL 


fi 


(c^)'sin“(0 )—0 5cos^(6 


+ sin‘(6 ) 


(61a) 


Similarly, solving (60b) for and taken into account (24) it follows that 


fisBL 

LBL 


(e^)‘cos^(ft)-0 5sin‘(8) 



(61b) 


Equations (61a, b) can be plotted in polar form in order to show the target's range (nor¬ 
malized to LBL) as it would be measured by a long and short baseline respectively 

D. RESULTS BASED ON SUM OF THE ABSOLUTES 

Following (55d) and (55f), the maximum normalized range LBL and SBL have been 
plotted in Figures 12, and 13 respectively, versus target orientation 6. T’nis was done for 
conditions a^=a^=l, c^=c^=500, and LBL=5 x SBL. For example, using 0=90° on Figure 
12 indicates that the maximum triangulation range is 250 x LBL. It is seen from Figure 13 
that the SBL maximum occurs at 0“ and it is 50 x LBL as expected. The combined effect 
is easily visualized in the cartesian representation shown on Figure 14. The solid line is the 
maximum triangulation range using the short baseline while the small circle plotting sym¬ 
bol is the corresponding long baseline range. As mentioned, the assumption inherent in this 
scheme is that the baseline producing the most precise estimate for range will automati¬ 
cally be selected with the use of "smart electronics." Therefore the larger of the short and 
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long baseline curves for represents the ideal capabilities for the dual baseline scheme 
This is shown on Figure 14 with a large circle used as a plotting symbol The more physi¬ 
cal polar representation for this solution is provided in Figure 15. In this figure can be seen 
the maximum range which is measured via the dual baseline scheme in any direction. 

E. RESULTS BASED ON QUADRATURE ADDITION OF TERMS 

Figures 16, 17, 18, and 19 are the counterparts for Figures 12, 13, 14, and 15 with 
the same data (i.e., a^=a^=\, c^==c^=500, and LBL=5 x SBL) but based on {61a), and 
(61b). It should be noted that these curves have a that as expected is approximately 
2' ■ higher (i.e., less conservative) compared with the previous case. 

F. OBSERVATIONS 

It is possible to obtain approximate analytic expressions for / LBL in terms of 
the target orientation. In the case of the more conservative sum of absolutes this can be 
done for arbitrary bearing precision ratios and (/. Figure 12, obtained using the sum of 
absolutes paradigm, shows that the maximum value for the long baseline range is approxi¬ 
mately 250 X LBL, where 250 is 50% of the long baseline sensitivity This can be com¬ 
pared with similar long baseline results for the quadrature addition paradigm shown on the 
Figure 16. Here the maximum value for the long baseline range is approximately 350 x 
LBL. where 350 is the 70% of the long baseline sensitivity. These observations, which 
show that the range scales with the sensitivity, can be generalized to convenient rules of 
thumb for calculating the maximum triangulated range for the long baseline system The 
details of this case are presented in Chapter VI. This suggest that for the condition c‘ 
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c c the general appearance of the two curves, independent of scale, will depend on the 
ratio of the long baseline to the short baseline. This issue is clarified also in greater depth 
in Chapter VI. 
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normalized range detected by long baseline (sum of absolutes) 
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ratio of the bearing errors in long baseline=1 normalized sensitivity of the long baseline =500 


Figure 12: Normalized range measured by the long baseline sensors for a^=l, and c’^=500 
(based on the sum of absolutes and obtained via program in Appendix E. D) 
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normalized range detected by the short baseline(sum of absolutes) 



ratio of bearing errors in short baselin=1 normalized sensitivity of the short baseline=500 


Figure 13: Normalized range measured by the short baseline for a^=l. c^=500, and 
SBL/LBL=0.2 (based on the sum of absolutes and obtained via program 
in Appendix E. E) 
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range detected by the dual baseline system in any direction(sum of absolutes) 



Figure 14: Cartesian plot of the normalized range measured by the short baseline, and the 
long baseline, showing also the range measured by the entire dual baseline sys-- 
stem, (based on the sum of absolutes and obtained via program in Appendix 
E. F) 
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range detected by the dual baseline system in any direction(sum of absolutes) 
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Figure 15 Plot of the nonmaJized range measured by the dual baseline system in any di¬ 
rection. (based on the sum of absolutes and obtained via program in Appen¬ 
dix E. G) 
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normalized range detected by long baseline(quadrature addition of terms) 
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Figure 16 Normalized range measured by the long baseline sensors for and c‘'=500 
(based on the quadrature addition & obtained via program in Appendix E H) 
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normalized range detected by the short baseline(quadrature addition of terms) 
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Figure 17: Normalized range measured by the shon baseline for a^=l, c’^=500, and 

SBL/LBL=0,2 (based on the quadrature addition and obtained via program 
in Appendix E. I) 
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Figure 18: Cartesian plot of the normalized range measured by the short baseline . and the 
long baseline, showing also the range measured by the entire dual baseline sy¬ 
stem (based on the quadrature addition and obtained via procram in Appendix 
E J) 
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range detected by the dual baseline system in any direction(quadrature addition) 
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Figure 19: Plot of the normalized range measured by the dual baseline system in any di¬ 
rection. (based on the quadrature addition and obtained via program in Appen¬ 
dix E. K) 
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V. PHYSICAL LIMITATIONS OF THE TRIANGULATION METHOD 


A. INTRODUCTION 

The polar plots of the normalized sensitivity that were presented in the previous sec¬ 
tions define limitations in the range measurement via the triangulation model. In this chap¬ 
ter the physical limitation due to the optical horizon is considered. 

B. GEOMETRIC AND OPTICAL HORIZON [REF. 18) 

The line at which the earth's curved surface apparently meets the sky is called the ho¬ 
rizon It can be observed very easily that the distance to the horizon increases as the height 
of an observer's eye increases above the earth surface. The dependence of sighting dis¬ 
tance on the height of the observer's eye is by itself an evidence that the surface of earth is 
curved As shown in Figure 20, two types of horizons can be distinguished: the geometric 
horizon that extends to the point H, and the optical horizon that extends to the point //'. 
The geometric horizon can be defined as the points where all the straight lines from a 
given observer tangentially graze the surface of the earth. However, due to refraction an 
observer can see a longer distance than that of the geometric horizon. This displaced hori¬ 
zon is the only one that can be sensed by an observer and it is called the optical horizon 
So, the refraction enables the observer to see over and beyond the earth's surface at H. 
The distance to which an observer can see depends on the height of the eye and the 
amount of refraction, and uncertainties in this latter quantity can be troublesome when 
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high accuracy is needed. The uncertainties arise from the fact that the line of sight to the 
horizon lies too close to the earth's surface. Because of differences in air and land or water 
temperature, abnormal temperature inversions or temperature gradients may occur in the 
lowest few meters of the atmosphere. Thus, the amount that the line of sight is bent is a 
variable and generally very difficult to predict, leading to a high level unpredictability in 
the optical horizon location. Therefore the more conservative estimation for the optical 
horizon, based on geometric.il consideration alone, will be followed. 

C. OPTICAL HORIZON GEOMETRIC ANALYSIS 

As may be understood from the Figure 20, no part of the earth's surface whose ele¬ 
vation angle is less than that of the optical horizon can be seen beyond this horizon How¬ 
ever, an elevated object such as an incoming missile, flying above the surface of the sea, 
can be seen beyond this horizon. The geometry for a such an object is shown on Figure 
21, and the foregoing analysis refers to an average steady-state condition of the air. From 
Figure 21, it is obvious that 


d = {Re + 1) • sin 6, + {Re + 5 ) • sin 0j 


(65) 


Again from the same figure 


COS0, 


/?,+/ 


(66a) 



(66b) 
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After (66), (65) becomes 



The geometric analysis ignores the effects of the refraction which can increase the distance 
to the optical horizon. Therefore (68) represents a conservative estimation for the optical 
horizon The relation (68) gives an approximation of the distance at which an object at 
height / above the surface of the earth can be seen by an observer or sensor at height i 
above the surface of the earth. Taking R^=6,37x10 * meters. (68) can be written as 

^(=3569.6(77 + 77) meters (69) 


where d, t. and s are expressed in meters. 
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D. EXAMPLES INCORPORATING THE GEOMETRIC HORIZON 


The dual baseline triangulation model permits a maximum triangulation range to be 
estimated within a maximum tolerable error criterion. In the case that the geometric hori¬ 
zon appears at a distance shoaer than the dual baseline prediction then the de-facto maxi¬ 
mum triangulation range is the geometric horizon. The combined limitations of the dual 
baseline triangulation system composed of the inherent and physical limitations can be 
modeled and represented with computer plots. Such plots are presented in Figures (22a, 
b) and (23a, b) and it can be seen that the limitations in the behavior of the system depends 
either on the distance from the baseline to the geometric horizon or the parameters of the 
system itself The data used to generate these plots are conveniently provided in Table 1 
It is noted that the specification </=</=! is consistent with the criterion that all bearing 
sensors have the same precision in bearing measurement. On the other hand the specifica¬ 
tion 6-^= c’‘'=c is consistent with the condition that the range measurement should be 
worked with the same precision limit regardless of the baseline used. Note that Figures 
22a and 23a are based on the sum of absolutes, while the Figures 22b and ?3b are based 
on the quadrature addition of terms On all figures, a^=a^=l, and c‘^=6'^=500. 

TABLE 1 DATA FOR THE FIGURES (22, 23) 


t(m) 


s (m) 


GH (m) 


LBL (m) 


SBL (m) 


Figures 22a, 23a 
Figures 22b, 23b 


20 

3 


10 22,480 

3 12,365 


100 

20 


* On all figures =0^=1 , and c^=c‘'’=500 
** Figures (22a), and (22b) are based on the sum of absolutes 
*** Figures (23a), and (23b) are based on the quadrature addition of terms 
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Figure 20 Horizon geometry and efTects of refraction (After Ref 18) 
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Figure 21 Geometrical horizon geometry in case that both sensor and target are above 
the surface of the earth 


64 






detected range is limited by the geometric horizon (sum of absolutes) 


90 

2.5e+004 



270 


height of the target above surtace^lOm height of the sensor above surf ace* 10m 


Figure Z2a. The maximum measured range is limited due co geometric horizon in 22A80 
m instead of 25.000m that can be detected by the system, (based on the sum 
of absolutes and obtained via program in .Appendix E. L) 



d@tected range is limited by the system (sum of absolutes) 

90 

1.4e+004 


range detected by the system 


geometric horizon 


180 



long baseline length-20m 
short baseline length*4m 240 

height of the target above surtace-3m 


270 


height of the sensor above surface* 


Figure Z2b The maximum deteaed range is limited to 5.000m due to the method, instead 
the I2.j65m which is the geometnc horizon (based on tne sum of absolutes 
and obtained via program in .Appendix E.L ) 


'3m 
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detected range is limited by the geometric horizon (quadrature addition) 


®°2.5e+004 



Figure 23 a The maximum measured range is limited due to geometric horizon in 22,480 
m instead of 35,000m that can be detected by the system, (based on the quad¬ 
rature addition of terms and obtained via program in .Appendix E, M) 
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detected range is limited by the system (quadrature addition) 


90 

1.4e+004 



height of target above surface*3m height of sensor above surface=3m 


Figure Zjb The maximum deteaed range is limited to 7,000m due to the method, instead 
the 12.j65m which is the geometric horizon (based on the quadrature addi¬ 
tion of terms and obtained via program in Appendix E.M) 
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VI. GENERALIZATION OF THE RESULTS FOR FIXED SENSITIVITY 


A. MOTIVATION 

As seen in the previous chapter the condition for fixed sensitivity leads to polar 
curves for maximum triangulation range using both a LBL and a SBL. The crossing point 
for these curves, 6 c.ross, defines the complimentary application domains for the respective 
baselines. In this chapter some simple guidelines for predicting Bcross- maximum SBL 
range, (at 6=0°), and the maximum LBL range, (at 6=90°), are derived. This is 
done for the quadrature addition of terms and the more conservative sum of absolutes 

B. A SIMPLE FORMULA FOR THE LBL & SBL CROSSOVER 
1. Analysis for the sum of the absolutes 

It is observed that in the crossover points of the LBL and SBL polar plots the 
normalized range in respect to LBL for both baselines has to be equal. That means 


Hlbl ^ Rgi 
LBL " LBL 


(70a) 


Taking into account (61), and (62), (70a) becomes 




(70b) 
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or 


sin‘(0) 




= COS' 


,|C-^)'+(l+fl^V 


( 1 + 0 *)- 


sm. 

LBL 


(70c) 


Furthermore (70c) becomes 


tan-(0) 



(71) 


Normally in (71), the terms of the form )“ are much bigger than the terms of the form 
(!+</''')• and dominate, so (71) can be written 

For the typical case when c^=c^=t'. Equation (72) becomes 

tan(e)=(^] (73a, b) 


and finally 


Across = Otan | ^ 


( 74 ) 


It is worth noting for the correction of (70) that if and c^=c^ then (74) is exact 
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2. Analysis based on the quadrature addition of terms 

Similarly to the previous case at the crossover points the normalized ranges of 
both baseline systems are equal. This implies that 



and finally after (63) and (65), (75) gives 


(c^) Mn^(9‘-)-0Scos^(9^) 


+ sin 


in‘(e^) = + sin-(e^) 


(76) 


and after (24), and (25), (76) become 


|f^rsin‘(9)-0 5cos^(9) . 2/n\ f (c^)*cos*(9)-0 5sm^(9) 

- 1 - + SU1 (0) = l - ^ -+ COS 




(77) 


or 


(c^^m^9W)5(l-5m=(9)] , f 5sm’(9) 

-;-+ sin‘(0) = 


^1 -sin‘(0) 


(Sr I (’8) 


After assuming that c‘^=c^=c and substituting sin^(e) with "y," in (78), it follows that 


0.5c*>'-0.25+ 1.25_>' = (0.5c'-0.5c‘>'+ 1-1. 



(79a) 


71 




and finally 


0.5c- + 1.25 + 0.5c-f^ 


+ 1.25 




5 + 0,5r| = 



(7%) 


Solving for "y," that is sin'(0), and after back substitution into (77), it follows that 


sin"(6) 


0 

l2MOSc=-K)5c=(^)'+125(§f)' 


(80) 


It can be observed at this point, that after dividing both terms of the ratio by c' and ne¬ 
glecting the very small terms that follow, (78) becomes 


sin‘(0) 




(81a) 


or 


fai) 

sin (e) = -p4^^ (81b) 


and finally 


r 


Across =osin 


ii(*7 


= atan 


SBl. 

LBL 


(82) 
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3. Observations 


After the analysis of both methods have been completed the very obvious result 
is that the crossover point depends only on the ratio of the two baselines and it is approxi¬ 
mately independent of the range sensitivity factor c. To demonstrate the above statement 
the more exact (80) and the approximate (82) crossover angle is plotted in Figure 24 for a 
fixed baselines ratio SBL/LBL=0.2 versus the variable normalized range sensitivity The 
approximate form (82) is nearly equivalent to (80). 

C. SIMPLE SCALING FORMULA FOR DUAL BASELINE MODEL RANGES 
1. Analysis for the sum of the absolutes 

As seen in Chapter IV the sensitivity-fixed triangulated ranges, normalized by 
the LBL, become maximum at the angle of 90“, measured with respect to each baseline. 
This maximum normalized range can be evaluated from expressions (61) and (62) for the 
long and short baselines respectively, as follow 


-JJJ—l0t=9(r - 


(1+a^)' 


(83a) 


and 


I _ I f SBL y 


After assuming that c^=c^=c, and c/=a^=l, (81) becomes 


(83b) 
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LfiL “ 



and 



(84a, b) 


2. Analysis for the quadrature addition of terms 

Following the same procedure described above (63), and (64) give 


^lAlunn) I 
LBL 



(85a) 


and 



Finally supposing c‘-=cf=c, (83) follows that 



(85b) 


(86a, b) 


As expected, a factor of 2''^ increase in the less conservative quadrature addition of the 
terms is seen by comparing with (84). 

3. Observations 

The most obvious result of the analysis above is that the scale of the 
curves increases linearly with the normalized range sensitivity. More specifically for the 
sum of absolutes paradigm the maximum ranges satisfy the following rules of thumb 
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= T xLfi/4m) and RsB!(m) = ^'x.SBL(m) 


(87a, b) 


Similarly, for the quadrature addition of terms paradigm the maximum ranges satisfy the 
following relations 


Ruii.{m) = -^xLBLim) and RsHL{m) = -^xSBL{m) (88a, b) 


For both paradigms the crossover angle is given by the rule of thumb 


Across = <?tan 



(89) 


The general characteristics of the dual baseline triangulated ranges and baseline crossing 
points are shown on Figure 25. 
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angle thela in degrees 


comparison between crossover angle approximations 


11.3103 


11.3102h- 


11.3102 


11.3101 


11.3101 


11.31 


11.31 


11.3099*- 

500 



650 700 750 800 850 

normalized range sensitivity 


Figure 24b: Crossover angle of long and short baseline based on the accurate expression 
(80) and the approximate expression (82) The discrepancy between the two 
expressions becomes obvious as the scale on the y a.xis goes to the third deci 
mal point, (via program in .Appendix E. N) 







case A 



(case A sum of absolutes , case B . quadrature addition of terms ) 


Figure 25 : Properties which characterize the performance of the dual baseline triangula¬ 
tion scheme The maximum normalized range which can be detected by the 
system is linearly dependent on the sensitivity. The crossover angle of the 
long and short baseline approximately depends only on the ratio of the two 
baselines 
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VII. CONCLUSIONS 


A. RESULTS OF THE WORK 

A dual baseline scheme based on the triangulation principle proposed for passive 
range measurements was presented and its performance was examined in this work In 
Chapter II the approximations on which the proposed model was based were examined, 
and the general characteristics of the triangulation principle were reviewed The conclu¬ 
sion was that the approximations R»BL, and R=R, =R 2 have not noticeably changed the 
accuracy of the model, and that for angles close to 0° and 180“ due to "geometric dilution" 
no range measurements can be performed. Furthermore it was shown that the precision in 
the range measurement was dependent on the precision in the bearing measurement. Lastly 
it was shown that the more distant the target to be triangulated is from the baseline the 
more inaccurate the measurement. In Chapter III, feasible modeling for a dual baseline tri¬ 
angulation scheme was examined. It was shown that two orthogonal baselines eliminate 
the "geometric dilution" problem and that, in general, the larger baseline will have a 
higher measurement accuracy. Finally, for further study of the characteristics of the dual 
baseline scheme two gauges were proposed. The first one is based on the quadrature addi¬ 
tion of the terms which express the range sensitivity of each sensor for each baseline, and 
the second one is approximate and is based on the sum of absolutes of the same terms It 
was shown in Appendix D that the latest one is more conservative with respect to the er¬ 
ror in the range measurement, and it can be used for arbitrary bearing precision of the 
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sensors. Instead, the less conservative quadrature addition of terms can be used only for 
the case when the sensors have the same bearing precision In Chapter IV the perform¬ 
ance of the dual baseline model was studied by the above mentioned gauges through the 
polar plots of the normalized range versus the target orientation angle. The conclusion 
from this chapter was that with the use of the conservative sum of absolutes the predicted 
maximum range is 2'' less than the range predicted by the quadrature addition of terms In 
Chapter V the physical limitation in the range measurements due to optical horizon was 
considered. The effect of refraction introduces a complicated situation in the prediction of 
the optical horizon. Thus, a more conservative estimation for the optical horizon based on 
the geometrical consideration alone was taken. It was shown in this chapter that the pre¬ 
dicted maximum range depends more generally on sensor and target height above the 
earth's surface, and the ratio of the two baseline lengths. Therefore the range measurement 
can be limited either due to the optical horizon or the dual baseline system In Chapter VI 
the general features of the polar plots that are characteristics of the dual baseline method 
were considered. It was shown that the maximum range scales with the range sensitivity of 
the baselines. In the case when the normalized sensitivity in both baselines is the same 
then the general appearance of the two curves will depend only on the ratio of the long 
baseline to the short baseline More specifically, the crossover angle of the two curves, to 
a high degree of approximation, depends only on the ratio of the two baselines, and it is 
independent of the range sensitivities. 
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B. FUTURE DEVELOPMENTS 


The dual baseline model was studied under the assumption that the mean bearing er¬ 
ror was zero, but in the real world that is not true. Thus future work is needed to enhance 
the model to handle systematic errors in the bearing measurements. This future work 
could employ the "circulation method" instead of triangulation. A second point for future 
development could be the implementation of the dual baseline model with some "smart 
electronics" which would switch between baselines at a performance crossover point , 
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APPENDIX A: LIST OF SYMBOLS 


A. SYMBOLS IN ENGLISH 

BL: the baseline (platform length) in the triangulation method (m) 
c. normalized range sensitivity in triangulation method 
normalized range sensitivity in long/short baseline 
d: distance to the geometrical horizon (m) 

LBL: the long baseline in the dual baseline triangulation scheme (m) 

R; the range from the middle of the baseline to the target (m) 

R^. radious of the earth (6.371 x 10® meters) 

Rlpl the range from the middle of the long baseline to the target (m) 

Rjgj^: the range from the middle of the short baseline to the target (m) 

R,,,: the range from the first second sensor to the target (m) 

SBL: the short baseline in the dual baseline triangulation scheme (m) 
s: height of the sensor above the surface of the earth (m) 
t: height of the target above the surface of the earth (m) 

B. SYMBOLS IN GREEK 

a : ratio of the mean values of the bearing errors in the triangulation method 
: ratio of the mean values of the bearing errors in long/short baseline 
A6„; : bearing error of the first/second sensor (rad) 
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6 : the target orientation angle measured from the middle of the baseline (degrees) 

6 ‘ ^ angle to the target measured from the middle of both baseline (degrees) 

6i; angle to the target measured from the first/second sensor (degrees) 

X| independent random variable having zero mean and unity standard deviation in 
long/short baseline in respect of first/second sensor 
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APPENDIX B: MATHEMATICAL PROOFS 


A. SINGLE BASELINE MODEL 
1. Single baseline model analysis 

The distance from the middle of the baseline to the target given by Equation 
(4) is repeated here 


R = ^ 


sindi 


Jl-(lysine,) 




1- frsme: 




(B. 1) 


The derivatives of (B.)) with respect 0,, and 0, are respectively 


■ 

/ 



“ 

cosHf 

1 

l-(fsine:] - 

sinBismBj 

-^(01.0;./?) 


v 



(B. 2) 


dR _ 



(] sinBisinfl; 


+ R(0,,B2.R) 


(B 3) 


where 


g(8i.e:,R) = 


( 




sinB: 



+ sin0 



(B 4) 


84 





and 


^ = + (B.5) 

The expressions (B.2), (B.3), and (B. 4) above are approximate because the derivatives of 
the small terms proportional to (BL/R)' on the RHS of Eq. (B 1) have not been included 
The terms 


can be factored out of Equations (B .2), and (B 4). Also the terms 


can be factored out between Equations (B 3), and (B.4), Then Equation (B 2), after Equa¬ 
tion (B 6), becomes 

f = (B.8a) 

where 

sinH; ) 

f\= ——(cosei+sin(Jicot(e:-ei)) (B 8b) 



sinSi ) 


sin(92-Bi) 


(B 7) 



Sin (62-61} 


sm02 j 


sin( 62 ”B|) 


(B 6) 
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and 


SmB; 


h 




sin(H2-B|) 




cot (62 -Bi)- 


sm Hi COS Hi 


(B. 8c) 


Similarly, Equation (B.3), after Equation (B 7), becomes 


dR 

dti 2 



(B 9a) 


where 



—0: - sin ejcot (9: - e,)) 


(B 9b) 


and 



sin 


001(02-0l) + 


V 


(Is ) 


(B, 9c) 


Equation (B. 5), after (B. 8), (B 9), and after dividing by (B. 1), becomes 


R 




^sinB:] (cosfli4«m9|Col(93-Bi)>+sin03 


f( 


fisuie,) 




COl( 02 - 9 l 


iin0i co«0 1 




r 



sinSj 

lysine, j +stn»i^ 



(ABiH 
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(cos 0 :-sinB 2 Coi(B:—0 ))-sin0 






COt(02—0j )-4— 




Sind I 


■Minfli Jl-(^s 


(Adz) (B 10) 


as cited in the main text. 

2. Single baseline model analysis (R»BL) 

A reasonable assumption that R»BL can be made, and it follows from Equa¬ 
tion (4)that 


R 


BL sinBj+sinSi 
2 sin(02-di) 


) 


The derivatives of (B. 11) with respect 0,, and 6, are 



BL I cos9i 
2 , sin ( 02 “®i) 


^(01.02)) 


and 


dR 


Bl I cos 9; 

2 I sm(9j-0,) 



(B. 11) 


(B 12) 


(B. 13) 
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where 


^(e,.0;) 




(B. 14) 


Equations (B. 12), and (B. 13) after (B. 14) can be re-expressed as follow 

^ = +C 0 t(e 2 -BiXsinBi +sine:)) (B. 15) 


and 


^-■f'i;;Hi^(“se:-cot(e;-ei)(sine:-sine:)) (B. 16) 

Combining Equations (B.15), and (B.16), and after dividing by (B.l 1) leads to 

f = (dS:5:+“'<9:-e')) w9.)+(=fS5r-cot(e:-e,)) (4e,) (b, h) 

as cited in the main text. 

B, DUAL BASELINE MODEL 

Following the substitutions suggested by Equations (23), (24), (26), and (27) the 
expressions for obtained for the various approximation cases of 

the dual baseline model analysis introduced in Chapter III. 
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1 . 


Dual baseline model analysis 


a. Long baseline (LBL) 


AR I 
R AH 


Ji-I^suiB; )* (cos8f+sinflfcoi(8{-8f ]]-«ine'^ 




sinSi j +smOf ^ 

'.-(^.0,8^) = 




(cosfl;+sin0icol(8;-e{ ) j-sinef^l-f^sinfl; ) coi(hJ-h{ 










+4m0 




(B 18) 


X:« 


6. baseline (SBL) 


^ J_ 
AH 


^ 1 ^ stn ef j ’ [cos0f+sin 0f cot (Sj-ef ) )+sin efjl- 

/ ...nSV 

t \ * i *10 07 to* S'? 

...r^s -.s'! « ' ' ' 

1 vt j 


smdi^ j 45in6 





xf+ 


I 


J ) (cos0f+5inflfcol| 

8f-0f))-sm0f^l-(^sm0f)' 

’ 

CO 

„s) 


sinBj^ 

r 

1- 

(^sm0f)’‘H,u,8fJl-(^sm0f) 



(B. 19) 


2. Dual baseline model (R»BL) 
a. Long baseline (LBL) 


(V=r 

\ ff AH j l, jme(-«in8j 


+cot ( 02 -et 



cos 0* 




■cot 





(B 20) 
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Short baseline (SBL) 


f cos 


V sin Hx-«in H‘| 


+ col B 




xf+ 


smBA+sinH', 


-cot B- -B, 


(B.2!) 
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APPENDIX C: EXACT EXPRESSION FOR THE RANGE 

A. GENERALITIES 

From Figure 2 the expression (3b), which relates the range of the target to the apex 
angles on the baseline, had been derived. Solving (3b) for R gives the exact expression 
that has to be used to solve the problem in the triangulation method. 

B. ANALYSIS 

In order to simplify the notation the following definitions are introduced 

ATI'S sin (02 - 01 ) 

^ s sin 02 
sw = sin 0 i 

Then from (3b) 



After substituting 

Z=(f) (C.3) 


(C. la) 
(C lb) 
(C. Ic) 
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symbolic processing using MAPLE [Ref. 19], produces 


4Z‘^.si ‘(-l +Z^sy^) = (Z'sy^ -IZ'^.sy^sM’^ +Z*.sit ‘ -.sn ')' (C. 4) 


and then again symbolic processing produces the two solutions 


2 j -Isx^sy^ -+ 4 J-sx'*.<i)’~su4(-] +sx^) 
2 -Asw’^sy^sx^ - ^ 3 ^"* + IsM'^sy^ - 5 ^“* 


(C 5a) 


- 1 -2sx^5V^ - ls\i4sx^ - 4 J-sx^sy^sH'^(-\ +5x^) 

(Z2)2 = i--^- 


2 -Asw^sy^sx^ — + 2^^^.^^ - w'* 


(C. 5b) 


It can be shown using MAXSYMA [Ref 20], that the root (C. 5a) is real and (C 5b) is 
imaginary. Therefore the exact expression for normalized range in terms of angles 0,. and 
0, is 



—45M'^^^sx^ —sy* + Isw'^sy'^ - 5M’'’ 
-Isx^sy^ — Isw^sx^ + 4 J-sx'^sy^sM'^i-l +sx^) 


(C.6) 
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APPENDIX D: RANGING ERROR IN TRIANGULATION METHOD 


A. GENERALITIES 

In single triangulation the target lies at the intersection of the two lines that run from 
the sensors. These two lines make angles with the baseline that vary due to mechanical 
limitations of the sensors. The variation in the angles, A0, and A6,. determine the varia- 
*'^n in the measured range of the target. 

B. APPROXIMATE RANGING ERROR ANALYSIS 

Figure 26 shows the geometry of range measurement in the single baseline triangula¬ 
tion scheme. The quadrangle KLMN is established due to variation in the angles. The lines 
DE, and FG have values AR, and AIL respectively, and approximately represent the un¬ 
certainties in the range of the target from the sensors. The line MK represents the total un¬ 
certainty due to the system and it has value AR. In first approximation the error 
quadrangle can be treated as a parallelogram and the following relation is derived from 
parallelogram on Figure 26 


ARq = ARparauelogram = ARp 


(D 1) 


or 


A/?,. = ^(Ay?,)- + (A/? 2 )- + 2|A^,| • |A/?;|cos(t) 


(D 2) 
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where "(f)" is the angle between AR, with AR^. 

For R»BL, as shown in Figure 27, <() = 0‘, and therefore 

cos(t) = 1 (D 3) 

Thus, it follows from Equation (D. 2) that 

Mp-M<lM,l + lM2l (D.4) 

and therefore via comparison with (D.2) it follows that 

V(A^i)' + (A/?2)‘ < |A/?, I +1 A/?; I (D.5) 

In fact, under the case that 


|AR,| = |AR:| 


(D.6) 


e g at a point of symmetry, and 


cos(t) = 1 


(D7) 


the sum of quadratures estimate will be 2 ''^ smaller than the sum of absolutes. So, in gen¬ 
eral, the sum of absolutes is a more conservative gauge for the lack of precision in the tri¬ 
angulation measurment. 
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BL 




Figure 26. Geometry for ranging error in single baseline triangulation 
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APPENDIX E: MATLAB CODES 


A. SINGLE BASELINE MODEL 

clear; 

b=input ('put the ratio RyBL='); 

y=l./b; 

a=input('put the ratio DTH2./DTH1-); 
step=pi./40; 
theta=0:step.pi; 
u=sin(theta); 
v=cos(theta); 
q=tan(theta); 

j=i /q; 

c=(b.^2+.5/^2+b.*v)/'.5; % c=R2/BL eq. 9a ; 

d=(b.^2+.5/'2-b.*v)/'.5; % d=Rl/BL eq. 9b ; 

thetal=asin(u *b./c); 
e=sin(thetal); 
f^cos(thetal); 
theta2=asin(u.*b./d); 
g=sin(theta2); 
h=cos(theta2); 
theta21=asin(l./d.*e); 
z=tan(theta21); 
k=l./z; 
l=(,5.*y), 
m=l.*g; 

n=(l-ni.^2).^.5; 
p=l,*e; 

q=(l-p.^2).-5; 
r=g.*q+e,*n; 
w=k-(1.^2.'»e,*f./q.^2); 
x=k+(l,''2.*g.*h./m.^2); 

B=n, *(f+e. *k)+g. *q. '^w; 

C=q *(h-g.*k)-e.'*n.*x; 

D=(f/(g+e)+k); 

E=(h/(g+e)-k); 

F=(sqrt(l-(y.*u).^2)./(y,*u)), 

% s represents the general solution for single baseline model eq 7 
s=((B./r).^2+((C./r),*a).^2).^.5; 


% 

eq. 10b; 

% 

eq. 10a; 

% 

eq 11 ; 


% 1=BL/2R; 

% m=(BL/2R)*sin(theta2); 
% p=(BL/2R)*sin(theta 1); 
% r=denom. of eq . 7 ; 
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% t represents the approximate solution (R»BL) eq 15 ; 
t=(D.^2+(E.*a).^2).^.5; 

% z represents the approximate solution (R»BL,R1=R2=R) eq 21 , 

z=((.5 *j+F)/'2+((.5.*j-F).*a).^2i/'.5; 

theta=theta.*180./pi; 

plot(theta,s,'*',theta,t,'o',theta,z,'+'); 

axis([0 180 0 .]); 

xlabel(Theta in degrees'); 
ylabelCjdelta R / (R * delta theta)('); 

titleCcomparison of single baseline model sensitivity curves'); 

gtext("^ for more accurate case (equation 7)'); 

gtext('o for R»BL approximation case (equation 15)'), 

gtextC + for R»BL & R=R1=R2 approximation case (equation 21)'); 

gtext('R/BL='); 

gtext('a='); 


B. DUAL BASELINE MODEL 


clear; 

rovlbl=input('put the ratio R/LBL-); % 

y=l./rovlbl; % 

rovsbl=input('put the ratio R/SBL-); % 

x=l./rovsbi, % 


a=input('put the ratio DTHL2/DTHL1-); 

aa=input('put the ratio DTHS2./DTHS1=); 

step=pi./40; 

theta=0;step.pi; 

u=sin(theta); 

v=cos(theta); 

m=( 1 -sin(theta)/'2)/'. 5, 


rovlbl=R/LBL; 

y=LBL/R; 

rovsbl=R/SBL; 

x=SBL/R; 


c=y *u; 


b=(l-c/'2)/'.5; 

d=(l-(-v.^2))/'‘;; 
g=x "^(-v); 
Hl-g/'2).^.5; 


k=v./u, % k=cot(theta); 

% c=(LBL/R)'''sin(theta); 

% b=( 1-((LBL/R)'*sin(theta)''2)^.5; 

% d=(J-(-cos(theta)^2)'',5, 

% g=(SBL/R)*(-cos(theta)); 

% Hl-((SBL/R)*(-cos(theta)))''2)^.5; 
% the following relation is for LBL model eq 33a ; 
l=((l/2 *k+b./c)/'2+((l/2.*k-b./c) *a)/'2)/'5, 

% the following relation is for SBL model eq 33b ; 

s=((l/2.*d/(-v)+f/g)/'2+((l/2 *d./(-v)-f./g),*aa)/'2)/'.5; 

theta=theta. * 180./pi; 

plot(theta,l,'o',theta,s,"*'); 

axis([0. 180. 0, 2.5*10^4]); 

titleC dual baseline model sensitivity curves'); 
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xlabeiCTheta in degrees'); 
ylabelCDelta R/(R * Delta Theta'); 
gtextCo for long baseline'), 
gtextC* for short baseline'), 
gtextCRTLBL-), 
gtext('R/SBL-); 
gtext('a='); 

gtext('ratio of bearing errors in long baseline='), 
gtext('ratio of bearing errors in short baseline^'), 

C. RANGE LIMITATIONS FOR FIXED TARGET ORIENTATION 

clear; 

maxrovbl=input(' input the maximum ratio of r/bl '); 
minrovbl=input(' input the minimum ratio of r/bl '), 
thetv=l:3; 

thetv( 1 )= input (' angle 1 in degrees other than 0,pi'); 
thetv( 1 )=pi/180*thetv( 1), 

thetv(2)= input (' angle2 in degrees other than 0,pi'); 
thetv(2)=pi/l 80*thetv(2); 

thetv(3)= input (' angle3 in degrees other than 0,pi'), 

thetv(3)=pi/180*thetv(3); 

nmax=input(' input the # points '); 

long=l nmax; 

longl = l nmax; 

long2=l :nmax; 

long3=l :nmax, 

longv=long'*thetv; 

difr= (maxrovbl-minrovbl)/(nmax-l); 

for k=l:3, 

thet=thetv(k); 

for i=l:nmax; 

rovbl(i)= minrovbl+difr*(i-1); 
rovbll=rovbl(i); 
v=sin(thet); 
z=tan(thet); 

x=l./z; 

long(i)= l/2.*x + (sqrt(l.((l /rovbll).*v)/'2))./((l Vrovbll) *v)); %eq 47 

longv(i,k)=long(i); 

end ; 

end; 

for i=l :nmax ; 
longl(i)=longv(i,l); 
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long2(i>="!.';ngv(i.2). 

lonp’'^-;-longv(i,3), 

fi' 

plot{rovbl,long 1 ,'o',rovbl,long2,'*’,rovbl,long3.'x') . 

title (' symmetrical analysis ') 

xlabel (' range/base-line ') : 

ylabel (' delta R/(R* delta theta) '), 

gtextCo for target approching at an angle -- degrees’), 

gtextC* for target approching at an angle — degrees'); 

gtextfx for target approching at an angle -- degrees'); 


D. LONG BASELINE POLAR PLOT (SUM OF ABSOLUTES) 

clear, 

%al=input('put the ratio of the bearing errors DTH2L/DTH1 L='), 
al=l; 

%cl=input('put the critical value of the normalized sensitivity^'); 
cl=500, 

theta=0 0 05:2*pi, 

% rlovlbl represents the eq.55d , 

rlovlbl=abs((sin(theta).^2 '^((cl/'2+(l+al)/'2.)/(l+al)/'2))/' 5); 
polar(theta, rlovlbl), 

titleCnormalized range detected by long baseline (sum of absolutes)'), 
gtextCratio of the bearing errors in long baseline='). 
gtext('normalized sensitivityof the long baseline=’); 
gtext('detected range'); 

E. SHORT BASELINE POLAR PLOT (SUM OF ABSOLUTES ) 

clear; 

%as=input('put the ratio of ratio of the bearing errors DTH2S/DTH1 S='); 
as=l, 

%f=input('ratio of the two baselines SBL(m)/LBL(m)='), 

1 ^ 0 . 2 , 

% cs=input('put the critical value of the normalized sensitivity^'); 
cs=500, 

theta=0:0,05;2*pi, 

% rsovlbl represents the eq 55f; 

rsovlbl=abs((cos(theta)/'2 *((cs/'2+(l+3s)/'2) /(l+as)/'2) 
polar(theta,rsovlbl); 

titleCnormalized range detected by the short baseline (sum of absolutes)'); 
gtext(’LBL=5 x SBL'), 

gtextfratio of the bearing errors in short baseline='). 
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gtext('normalized sensitivity of the short baseline-), 
gtext('detected range'); 


F. LBL & SBL CARTESIAN PLOT (SUM OF ABSOLUTES) 

clear, 

%cl=input('critical value of the normalized sensitivity in long baseline='), 
cl=500, 

%al=input('ratio of the bearing errors in LBL (DTH2/DTH1)='); 
al=l; 

%f==input('ratio (SBL(m) /LBL(m)='); 
f^0,2; 

%cs=input('critical value of the normalized sensitivity in short baseline ='); 
cs=500; 

%as=input('ratio of the bearing errors in SBL (DTH2/DTH1 )='), 
as=l, 

theta=0:0.05:2*pi; 

rlovlbl=abs((sin(theta)/'2.*((cl.^2+(l+al).^2)./(l+al)/'2))/'.5); % eq 55d ; 

rsovlbl=abs((cos(theta)/'2 *((cs. ^2+( 1 +as)/'2)./( 1 +as)/'2). *.. 
f''2)/^.5); % eq. 55f, 

rmax=max(rlovlbl,rsovlbl); 

%polar(theta,rmax), 
theta=theta.’* 180 /pi; 

plot(theta,rlovlbl,'.',theta,rsovlbl,'-.',theta,rmax,'o'); 

%gtext('al=2'); 

%gtext('cl=500'); 

%gtext('SBL/LBL=0.2'); 
xlabeK'theta in degrees'); 

%ylabel('rlbl/lbl, rsbl/lbl, rmax/lbl'); 

title(’cartesian plot of the range detected by the dual baseline system (sum of 
absolutes)'); 

gtext('ratio of the bearing errors in both baselines^l'); 
gtext('normalized sensitivity in both baselines=500'); 
gtext('LBL=5 x SBL'), 
gtext('range detected by short baseline'); 
gtext('range detected by long baseline'); 
gtext('range detected by the system'); 

G. LBL & SBL POLAR PLOT (SUM OF ABSOLUTES) 
clear, 

%cl=input('critical value of normalized sensitivity in long baseline='); 
cl=500. 
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%a!=input('ratio of the bearing errors in long baseline (DTH2L/DTH1L)='); 
al=l; 

%f=input('ratio (SBL(m) /LBL(ni)-); 
f=0.2, 

%cs=input('critical value of normalized sensitivity in short baseline ='), 
cs=500, 

%as=input('ratio of the bearing errors in short baseline (DTH2S/DTH1S)='); 
as=l; 

theta=0;0.05:2*pi; 

rlovlbi=abs((sin(theta)/'2.*((cl.^2+(l+al)/'2)./(l+al)/'2))/'.5); % eq. 55d ; 

rsovlbl=abs((cos(theta)/'2. *((cs/"2+( l+as)/'2)./( 1 +as).'^2). . 

f^2).^.5); % eq. 55f; 

nnax=max(rlovlbl,rsovlbl); 

polar(theta,rlovlbl,'.theta,rsovlbl,'- ',theta,rmax,'o'); 

title('range detected by the dual baseline system in any direction(sum of absolutes)'), 
gtext('ratio of bearing errors in both baselines= 1'); 
gtext('normalized sensitivity in both baselines=500'); 
gtext('LBL-5 x SBL'); 

H. LONG BASELINE POLAR PLOT (QUADRATURE ADDITION OF TERMS) 

clear; 

theta=0:0.05:2*pi; 

%cll=input('critical value of the normalized sensitivity in long baseline='); 
cll=500; 

% rlovlbl represents eq. 61a ; 

rlovlbl=abs(((l/2.*(cll.^2,*sin(theta).^2-0.5.*cos(theta).^2)+sin(theta).^2)).^.5); 
polar(theta, rlovlbl); 

titleCnormalized range detected by long baseline (quadrature addition of terms)'); 
gtext('ratio of the bearing errors in long baseline=r); 
gtext('normalized sensitivity of the long baseline=500'); 
title('detected range'); 

I. SBL POLAR PLOT (QUADRATURE ADDITION OF TERMS) 
clear; 

theta=0:0.05:2’'‘pi; 

%css=input('critical value of the normalized sensitivity in short baseline='); 
css=500; 

%f=input('ratio SBL(m) /LBL(m)-); 
f=0.2; 

% rsovlbl represents eq 61b, 

rsovlbl=abs(((l/2.’*(css/'2.*cos(theta)/'2-0.5 *sin(theta)/'2)+ 
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cos(theta).'^2)/‘f/'2).^0.5); 

polar(theta,rsovlbl); 

gtext('ratio of the bearing errors in short baseline=r); 
gtext('normalized sensitivity of the short baseline=500'); 

title('normalized range detected by the short baseline (quadrature addition of terms'), 
gtext('LBL=5 x SBL'); 
gtext('detected range'); 

J. LBL & SBL CARTESIAN PLOT (QUADRATURE ADDITION OF TERMS) 

clear; 

%cll=input('critical value of the normalized sensitivity in long baseline='); 
cll=500; 

%ff=input('ratio SBL(m) / LBL(m)-); 
fF=0.* 

%css=input('critical value ot the normalized sensitivity in short baseline='); 
css=500; 

theta=0;0.05:2*pi; 

% rlovlbl represents eq. 61a ; 
rlovlbl=abs(((l/2 *(cll/'2.*sin(theta)/'2-... 

0.5,*cos^^theta)/'2)+sin(theta).^2))/O.5); 

%rsovlbl represents eq. 61b; 

rsovlbl=abs((( 1 /2. * (css. ^2. * cos(theta). ^2-0.5. * sin(theta). ^2)+,.. 

cos(theta).''2).*fF.^2).^0,5); 

rmax=max(rlovlbl,rsovlbl); 

plot(theta,rlovlbl,'.',theta,rsovlbl,'-.',theta,rmax,'o'); 
xlabel('theta in degrees'); 
ylabel('rlbl/lb!, rsbl/lbl, rmax/lbV); 

title('cartesian plot of the range detected by the dual baseline systerr.(quadrature 
addition)'); 

gtext('ratio of bearing errors in both baselines =1'); 
gtext('normalized sensitivity in both baselines=500'); 
gtext('LBL=5 x SBL'); 
gtextC range detected by short baseline'); 
gtext('range detected by long baseline'); 

K. LBL & SBL POLAR PLOT (QUADRATURE ADDITION OF TERMS) 
clear; 

%cll=input('critical value of the normalized sensitivity in long baseline='); 
cll=500; 

%fF=input('ratio SBL(m) / LBL(m)-), 
fF=0.2; 
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%css=input('critical value of the normalized sensitivity in short baseline='); 
css=500; 

theta=0:0.05:2*pi; 

% rlovlbl represents eq. 61a ; 

rlovlbl=abs((( l/2.*(cll.^2.*sin(theta).^2-0.5 *cos(theta)/'2)+sin(theta)/'2))/^0.5). 

% rsovlbl represents eq. 61b ; 

rsovlbl=abs((( 1 /2. *(css, ^^2, *cos(theta).^2-0.5. *sin(theta). ^2)+... 

cos(theta)/^2).*fr.^2).^0.5); 

rmax=max(rlovlbl,rsovlbl); 

polar(theta,rmax); 

theta=theta. * 180./pi; 

titleCrange detected by the dual baseline system in any direction (quadrature 
addition)'); 

gtext('ratio of bearing errors in both baselines=]'); 
gtext('normalized sensitivity in both baselines=500'); 
gtext('LBL=5 x SBL'); 
gtext(’ass=r); 

gtext('range detected by the system'); 

LIMITATIONS IN THE MEASURMENT OF THE RANGE (S. A) 

clear; 

% S. A is for the sum of absoluets; 

%al=input('ratio of the bearing errors in the long baseline( DTH2L/DTH1L)='); 
al=l; 

%cl=input('critical value of the normalized sensitivity in long baseline='); 
cl=500; 

%lbl=input('length of the long baseline(m)-); 

lbl= 100/20; 

theta=0:0.05;2*pi; 

% rl represents the eq. 55d; 

rl=abs((sin(theta).^2 *lbl.''2.*((cl.^2+(l+al).''2)./... 

()+al).^2)).''0.5); 

%as=input('ratio of the bearing errors in short baseline (DTH2S/DTH1S)-); 
as=l; 

%cs=input('critical value of the normalized sensitivity in short baseline='); 
cs=500; 

%sbl=input('length of SBL(m)='); 
sbl=20/4; 

% rs represents the eq. 55f; 

rs=abs((cos(theta).''2. ’*sbl.''2. *((cs.''2+( 1 +as).^2)./... 

(l+as)^2)).^.5), 

rmax=max(rl,rs); 
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% re represents the radius of the earth (6.371 x 10^6 meters); 
re=6.371* 10.^6; 

%s=input('height of the sensor above surface (m)-); 
s=10/3; 

%t=input('height of the target above surface (m)-); 
t=10/3; 

% d represents the eq. 68; 
d=(2.*re.*t).'^.5+(2.*re.*s).'^.5; 

% r represents the geometrical horizon; 
r=(d. ^2. *cos(theta). ^2+d. ^2. * sin(theta). ^2). 5; 
if d>rmax 
polar(theta,r); 
hold on 

polar(theta,rmax); 
else 

polar(theta,rmax); 
hold on 
polar(theta,r); 
end; 

title('detected range is limited by the geometric horizon / by the system (sum of 
absolutes)'); 

gtext('long baseline length® 100/20m'); 
gtext('short baseline length =20/4 m'); 
gtext('height of the target above surface=I0/3m’); 
gtext('height of the sensor above surface =l0/3m'); 
gtextCgeometric horizon'); 
gtext('range detected by the system'); 
gtext('detected range'); 

M. LIMITATIONS IN THE MEASURMENT OF THE RANGE (Q. A) 

clear; 

% Q.A is for the quadrature addition of terms; 

%cll=input('critical value of the normalized sensitivity in long baseline®'); 
cll=500; 

%lbl=input('length of the long baseline(m)='); 
lbl= 100/20; 

%sbl=input('length of the short baseline(m)-); 
sbl=20/4; 

%css=input('critical value of the normalized sensitivity in short baseline®'); 
css=500; 

theta=0:0.05;2’'‘pi; 

% rl represents eq. 61a; 
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rl=abs((( 1/2. *(cIl/'2, *sin(theta)/'2-0.5. *cos(theta)/'2)+... 
sin(theta).^2).*lbl.^2).^0,5); 

% rs represents eq. 61b; 

rs=abs(((l/2.*(css/'2.*cos(theta)/'2-0.5.*sin(theta)/'2)+... 

cos(theta).^2).*sbl.^2)/'0.5); 

rmax=max(rl,rs); 

% re represents the radious of the earth (6.37 lx 10^6 meters); 
re=6,371* 10,^6; 

%s=input('height of sensor above surface (m)='); 
s=10/3; 

%t=input('height of terget above surface (m)='); 
t=10/3; 

% d represents eq. 68; 
d=(2.*re.*t).^.5+(2.*re.*s).^.5; 

% r represents the geometrical horizon, 
r=(d. ^^2. * cos( theta). '^2+d. ^2. * sin(theta). ^2). 5; 
if d>rmax 
polar(theta,r); 
hold on 

polar( theta, rmax); 
else 

polar(theta,rmax); 
hold on 
polar(theta,r); 
end; 

titleCdetected range is limited by the geometric horizon/by the system (quadrature 
addition)'); 

gtext('long baseline length=100/20m'); 
gtext('short baseline length=20/4m'); 
gtext('height of the target above the surface = 10/3m'); 
gtext('heigt of the sensor above the surface= 10/3m'); 
gtextCgeometric horizon'); 
gtext(’range detected by the system'); 
gtext('detected range'); 

N. CROSSOVER ANGLE 
clear; 

% f=input('ratio of the two baselines SBL(m)/LBL(m)='); 
f=0.2; 

% c is the normalized sensitivity of both baselines which assumed to be the same; 

% clower= input ('lower limit of the normalized sensitivity='); 
clower=500; 
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% cupper+input ('upper limit of the normalized sensitivity ='); 
cupper=1000; 

%cstep=input ('step of normalized sensitivity'); 
cstep= 10; 

c=clower: cstep xupper; 

thetal=atan(f*(c./c)); % eq. 82 ; 

theta2=asin(sqrt((0.25+0.5. *c/^2. *f "'2+. 

f^2),/(1.25+0.5.*c.'^2+0.5.*c.^2.*f^2+f^2))); % eq. 80 ; 

theta 1 =theta 1. * 180/pi; 

theta2=theta2. * 180/pi; 

plot(c,thetal); 

hold on; 

plot(c,theta2,'.'); 
axis([500 1000 11.3 11.4]); 
xlabel('normalized range sensitivity'); 
ylabel('angle theta in degrees'); 
title('comparison crossover angle approximations'); 
gtextCaccurate cross angle (eq. 80)'); 
gtext('approximate cross angle (eq. 82)'); 
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